Discrete Wave Turbulence 
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In this Letter we present discrete wave turbulence (DWT) as a counterpart of classical statistical 
wave turbulence (SWT). DWT is characterized by resonance clustering, not by the size of clusters, 
i.e. it includes, but is not reduced to, the study of low-dimensional systems. Clusters with integrable 
and chaotic dynamics co-exist in different sub-spaces of the k-space. NR-diagrams are introduced, 
a handy graphical presentation of an arbitrary resonance cluster allowing to reconstruct uniquely 
dynamical system describing the cluster. DWT is shown to be a novel research field in nonlinear 
science, with its own methods, achievements and application areas. 



1. Introduction. In 1960-1980th a great volume of 
research has been performed in the area of wave inter- 
actions [H, comprehensive account of theory and exper- 
iment can be found in Q. The notion of wave kinetic 
equation was introduced using some statistical as- 
sumptions, general methods for deriving kinetic equa- 
tions and their stationary solutions (energy spectra) were 
developed and statistical wave turbulence (SWT) theory 
was founded [4], with finite-size effects left aside. Their 
preliminary studies were performed in [1, 0], where it 
was established that nonlinear resonances are divided 
into dynamically independent, non-intersecting clusters. 
Explicit constructing of resonance clustering became the 
immediate challenge of great intricacy because no analyt- 
ical methods for solving resonance conditions were known 
(the problem in its general form is equivalent to Hilbert's 
10th problem [tJ); brute- force computer computations do 
not help either, while integers to be dealt with, are too 
big. The problem has been recently solved For com- 
plete resonance set, classical resonance curves ^ are not 
anymore a suitable representation, of course. Instead, 
representation by a hyper-graph on a plane was intro- 
duced in [10] (for 3- wave systems), which allows to ex- 
tract uniquely the dynamical system describing each clus- 
ter. The study of resonant clustering yielded a model of 
laminated turbulence [11] which gives kinematic explana- 
tion of co-existence of STW and DWT in wave turbulent 
systems. These and other results, both for 3- and 4- wave 
systems, are reviewed below, as well as their physical rel- 
evance and application areas. 

2. Why are predictions of SWT theory often not 
corroborated? The SWT theory assumes weak nonlin- 
earity, randomness of phases and infinite-box limit, i.e. 
the resonance broadening Q is greater than the spacing 
5 U between adjacent wave modes 

n > (<90/<9/c)2tt/L, (1) 

where L is the box size. Additionally, existence of in- 
ertial interval (k Q ,ki) is assumed, where energy input 
and dissipation are balanced. If k ^> k±, dissipation sup- 
presses nonlinear dynamics (see FigJTJ the upper panel). 
If k <C &o, finite-size effects take place which are due to 
boundary conditions and should be regarded separately. 
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FIG. 1: Schematic representation of classical SWT (upper 
panel) and both SWT and DWT (lower panel). Here E, k and 
(ko,ki) are notations for energy, wave- number and inertial 
interval correspondingly. 



In the last decade, this standard view proved to be 
incomplete while finite-size effects are well observable 
within the inertial interval ([HI, EH, 03 etc.). As it 
was shown in [ll|, this is due to the notorious small 
divisor problem: each small divisor leaves a gap in the 
power spectrum, shown as an empty circle in FigUl lower 
panel. The radius 1Z of each circle can be computed us- 
ing the Thue-Siegel-Roth theorem [l9[ for a big class of 
dispersion functions, which covers in particularly vari- 
ous types of water waves, oceanic planetary waves, drift 
waves in laboratory plasma etc. Gaps corresponding to 
exact and quasi-resonances are shown by yellow squares. 
Some gaps, shown by blue diamonds, correspond to dis- 
crete modes that do not take part in resonances and just 
keep their initial energy on the corresponding time-scale 
[5[. Accordingly, 3 types of solutions of ([2]),(j3]) must be 
distinguished: 

I: exact resonances, Q = 0; 

II: quasi-resonances, < 1Z < Q<(d0 / dk)2n / L; 
III: approximate interactions, Q ^> (d0 / 'dk)2n / L . 
Types I and II are the subject of DWT theory while Type 



FIG. 2: Examples of clusters in a 3-wave system. Upper panel: Topological structures of each cluster, Lower panel: 
NR-diagrams of the clusters shown in the upper panel. 
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FIG. 3: Examples of clusters in 4-wave systems. Left: Topological representation of a quartet. Middle: NR-diagram of a 
V-scale-scale cluster. Right: NR-diagram of a mixed cascade cluster, with one E- angle- scale, one VD-scale-scale and two 
V-scale-angle connections. 



Ill is covered by SWT theory. Resonance clusters of the 
first two types include low-dimensional systems but are 
not reduced to them: some clusters consist of hundreds 
or even thousands connected triads or quartets. These 
considerations justify the two-layer model of laminated 
turbulence, consisting of discrete and continuous (statis- 
tical) layers, presented in where the notion of DWT 
was first introduced. 

3. Exact resonances. 3-wave and 4-wave resonance 
conditions have the form 

01 + 02 = 03 + ^, ki+k 2 =k 3 , and (2) 
01 + 02 = 03 + 04 + ^, ki+k 2 =k 3 + k 4 . (3) 

correspondingly. Minimal possible clusters in three- and 
four-wave systems are triads and quartets correspond- 
ingly, they are called primary clusters. Cluster of arbi- 
trary structure can be constructed from primary clusters, 
as well as its dynamical system. All constructions below 
for concreteness are presented for 2D- wave vectors. 
3.1. Resonance clustering in three-wave systems. 
To construct complete resonance clustering in a given 
spectral place, three consequent steps have to be per- 
formed. 

At the first step, each resonant mode with wave vector 
k = (m, n) is presented by a node of the two-dimensional 
integer lattice (m, n) and every three nodes constituting 
a solution of ([2]) are connected by lines. The resulting 
graph is called geometrical structure (see [loj, Fig.l, up- 
per panel). 

At the second step, topological structure is extracted, 
which consists of all connected components found in the 



geometrical structure. Primary clusters are shown as tri- 
angles and their dynamical systems have standard form 

B\ = ZB^Bs, B2 = ZB^Bs, B3 = —ZB1B2, (4) 
differing only in magnitudes of the coupling coefficient 
Z. Due to Hasselmann's criterion of nonlinear wave in- 
stability [15], dynamical properties of the 03-mode with 
the highest frequency differ from those with smaller fre- 
quencies. We follow [16[ denoting 03-mode as A- and two 
other as P-modes. 

At the third step, a hyper- graph with marked arcs, is 
constructed whose (hyper-)vertexes are triangles con- 
structed at the previous step (see Fig|2j upper panel). 
To simplify graphical presentation for bigger clusters, we 
prefer NR-diagram presentation (NR for nonlinear res- 
onance), without letters and arrows but with two types 
of half-edges instead: bold for A-mode and dotted - for 
P-mode. 

Given a NR-diagram, the corresponding dynamical 
system is constructed by coupling 4, 3 or 2 systems of 
the form (|4|) and equaling appropriate Bi and Bj. For 
instance, the dynamical system of the two-triad cluster 
shown in Figj2] (it is unique, up to the change of indices 
1 <-» 2 and 4^5) reads 

B\ = ZB2B3, B2 = ZB\B%, Bs = —ZB1B2, 
< B4 = ZB^Bq, £? 5 = ZB%B§, Bq = —ZB4B5, => 

B\ = B4 (PP-connection), B3 = B$ (AP-connection) 
J B\ = ZB^Bs + ZB^Bq, B2 = ZB1B3, 

\Bs — —ZB1B2 + ZB^Bq, Bq = —ZB1B3, 
We have studied resonance clustering in data sets com- 
puted for various 2D- wave systems (different types of wa- 
ter waves, oceanic and atmospheric planetary waves, drift 



3 



waves, etc.) in the model spectral domain |m|, \n\ < 10 3 . 
In average, about 80 — 90% of all clusters consist of only 
one or two connected triads, though bigger clusters of a 
few thousand connected triads were also observed. Some 
clusters have integrable dynamics 17, 0], depending 



on coupling coefficients Z and/or on initial conditions 
and/or connection type, while dynamics of other clusters 
is chaotic. 

A fact of the major importance is: chaotic behavior is 
observable already in the systems consisting of two triads 
(AA-connection via one vertex, Z/Z = 0.75, work in 
progress). On the other hand, clusters of special form 
be integrable for arbitrary finite number of triads 
This means that clusters with chaotic and integrable 
dynamics co-exist, in different sub- spaces of the k-space. 

Summarizing, an arbitrary cluster appearing in a 3- 
wave systems can be represented as a NR-diagram with 
1 type of vertices and two types of marked edges: a sin- 
gle line for vertex-connection and a double line for edge- 
connection. Each edge is marked as A A, AP or PP. 
The presentation is sufficient to reconstruct uniquely the 
dynamical system of an arbitrary cluster. 
3.2. Resonance clustering in four- wave systems. 
Any 3-wave resonance generates energy transport over 
the scales in the k-space; this is not the case for 4-wave 
systems where three types of energy transport are ob- 
served: a) over scales, b) over angles, and c) mixed cas- 
cades, including both scale- and angle-resonances [T^ . 
Correspondingly, our NR-diagram representation of a 4- 
wave system has two types of vertices - circles for angle- 
resonances and squares for scale-resonances. Another 
important dynamical difference between 3- and 4-wave 
systems is that for quartets, no generic criterion of in- 
stability is known [15[; however, for the special case of 
0i = 02 or 03 = 04, the 3-wave criterion can be used, 
of course. To keep information about the vertices corre- 
sponding to the terms on the left and right side of ([3]), in 
geometrical representation each quartet consists of two 
edges (bold lines) connecting vertices from the same side 
of ([3]) and two diagonals (dashed lines). Accordingly, in 
NR-diagram representation, connections are possible via 
a vertex (a single bold line), edge (double bold), edge- 
diagonal (double, bold and dashed), diagonal-diagonal 
(double, dashed and dashed), notated as V-, E-, ED- 
and DD-connections correspondingly (see Fig|3|). Notice 
that connection types A and P in a 3-wave system, shown 
in corresponding NR-diagrams, define dynamical charac- 
teristics of a connection. In a 4-wave system, connection 
types V, E and D do not hold dynamical information in 
the general case, only in some particular cases. 
4. Diagrams for DWT versus diagrams for SWT. 
Various diagram techniques, beginning with Feynman di- 
agrams, are widely used for description of STW. NR- 
diagram is a tool to represent DWT; main differences 
between these two types of diagrams can be formulated 
as follows: for a resonance cluster look somewhat similar 



to the Feynman diagrams known in quantum mechanics. 
The differences are as follows: 1) In a Feynman diagram 
each vertex represent a particle, which corresponds to a 
single wave in the topological representation of a clus- 
ter, while in NR-diagram vertexes are primary clusters 
- resonant triads or quartets, depending on the order of 
resonance in the wave system under consideration. 2) 
Only a sum of all Feynman diagrams represents possi- 
ble interactions of a given particle with other particles. 
On the contrary, one NR-diagram represents complete 
resonance clusters. 3) A Feynman diagram does not al- 
low computing the amplitudes of the scattering process, 
it only gives a contribution corresponding to one term 
in the perturbation expansion. On the contrary, NR- 
diagram allows to reconstruct uniquely dynamical system 
whose solutions are the amplitudes of resonance cluster. 
5. Energy transport in DWT. The two-layer model 
of laminated turbulence gives only kinematic explanation 
for the generic co-existence of SWT and DWT. Numer- 
ical simulations with Hamiltonian dynamical equations 
[lH demonstrate the same. But we do not know any- 
thing about dynamical mechanisms underlying this co- 
existence. The problems to coop with are the follow- 
ing. In SWT, energy flow is described by energy spectra 
k _a which means, in particular that a) waves with wave- 
lengths of the same order also have energies of the same 
order; b) the magnitude of the initial energy of the sys- 
tem, Eo , is not important - only the fact that Eq = const 
is used; and c) wave phases are random. Statements a)- 
c) do not hold in DWT. Indeed, in many 3-wave systems 
waves with wave-lengths of order k and k 2 can interact 
resonantly ("relaxed" locality, 0), in 4-wave systems - 
waves with arbitrary big difference in wave-lengths can 
resonate (no locality, [8[). This means that the magni- 
tudes of modes' energies are not defined by the scales in 
THE k-space but by the initial values of modes ' energies 
and phases. Importance of phases in low-dimensional 
systems has been realized long ago [21[ but no analytical 
expression for phases was known. The most compact ex- 
pressions for energy and phase evolution in DWT regimes 
known to us presently read: 



EmodeiT^j 



cot I <p{T) | 
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where Em d e is THE energy of a resonant mode in a 3- 
wave system and (p is the dynamical phase [l7| which is 
a combination of initial phases corresponding to the res- 
onance conditions ([2]). Here Z is the coupling coefficient 
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from (|3J) and depends only on wave- numbers; /i, r, to and 
q are known expressions, including elliptic integrals, el- 
liptic functions, their nomes and modula. They depend 
on wave- numbers, initial energy and initial energy dis- 
tribution within a triad: not only E$ is important but 
also the part of it contained in each mode. Expressions 
similar to ([5]) , (|6]) can probably be obtained for a 4- wave 
system, along the lines given in [22[ • An important fact is 
that for a fixed resonant triad and fixed initial conditions, 
Z, /i,r, to,q are constants, i.e. ([5j), (|6]) can be used di- 
rectly for computing exactly energy and phase evolution 
of modes in primary clusters. These formulas can also be 
used for computing approximately energy and phase evo- 
lution in clusters of more general structure (for further 
discussion see [23], Chapter 4). 
6. Summary. 

• O. Phillips pointed out in [24J that SWT theory has 
reached its limitations and "new physics, new mathemat- 
ics and new intuition is required" for understanding dis- 
crete effects in wave turbulent systems. DWT theory 
gives reliable description of these effects and is indeed a 
natural counterpart to the classical SWT. In particular, 
this means that the classical notion of wave interactions 
is not sufficient and might even be misleading for de- 
scribing a generic wave turbulent system while it does 
not allow to distinguish between discrete and statistical 
regimes. Notions of exact resonances, quasi-resonances 
and approximate interactions should be used. 

• DWT is characterized by the clustering itself, and 
not by the number of modes in particular clusters which 
can be fairly big. This means that though chaotic dynam- 
ics can be accounted for in DWT, it can not be generally 
described by kinetic equations: each cluster has its own 
independent dynamics, sometimes integrable, sometimes 
not. Clusters of a given form generate persisting patterns 
observable in laboratory experiments ( [25[ , gravity water 
waves) and in measured data (0, atmospheric plane- 
tary waves). Their existence might give a better expla- 
nation of such well-known phenomena as Benjamin- Fair 
instability [ItJ (for a profound discussion, see [1^]) or 
freak waves (work in progress). 

• Numerical simulations [29[ show (for gravity water 
waves) that interactions of Types I and II are observable 
on the linear time scale, not on the time scale 0(t/e 4 ) as 
predicted by SWT theory. This indicates that DWT the- 
ory might give more appropriate foundation for a short- 
term forecast. 

• In some physical systems, only SWT is observable 
[30( ] (on-cite measurements of tidal currents), in other - 
only DWT [T3] (laboratory experiments). Co-existence 
of SWT and DWT has been demonstrated in numerical 
experiments [l4| (gravity water waves) where this regime 
was called mesoscopic wave turbulence. Including ad- 
ditional physical parameters can yield transition from 
SWT to DWT [3l| (capillary water waves, with and with- 
out rotation). Fluid mechanics examples do not exhaust, 



of course, the application areas of DWT, which include 
(similar to classical SWT) biology, medicine, astronomy, 
chemistry, sociology, etc. The problem of the utmost im- 
portance presently is to gain more understanding about 
SWT «=> DWT energy transport. 
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